Pinning of an Abrikosov vortex on a small cylindrical cavity: A Ginzburg-Landau 

approach 
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Within the Ginzburg-Landau theory we consider Abrikosov vortex pinning on a columnar defect 
with the characteristic size of the cross-section D much smaller than the coherence length f. We 
present an extension of the electrostatic analogies method, which proved to be useful for calculations 
of the pinning force for large cavities (D 3> £), to the case of small defects (D <C £). The pinning 
potential for an elliptic cavity is derived analytically. Also, we determine the depinning current for 
a circular defect. 

PACS numbers: 74.25.Wx, 74.20.De 



I. INTRODUCTION 

The possibilities for practical applications of type-II 
superconductors depend crucially on the quality of vor- 
tex pinning structures, which can be embedded into these 
materials. Columnar defects proved to be the most effi- 
cient pinning centers. Nowadays, various techniques exist 
allowing to create disordered arrays of such defects^ as 
well as regular defect lattices^ 

The first theoretical study of vortex pinning on a cylin- 
drical defect has been carried out by Mkrtchyan and 
Schmidt.— In their paper the London equation has been 
solved exactly for a vortex interacting with a cavity in 
the form of a circular cylinder. The pinning force has 
been analyzed in detail for the cavity radius a satisfiyng 
the condition («a<A, where £ is the coherence length 
and A is the London length. Later ^ this analysis has been 
extended to the case of large cavities with a ^> A. Buzdin 
and Feinberg 5 pointed out that London screening can be 
neglected in a large range of fields in extreme type-II 
superconductors. This observation allowed them to es- 
tablish an electrostatic analogy and to simplify consider- 
ably the solution for a vortex interacting with a circular 
cavity: it has been demonstrated that the full magnetic 
field can be presented as the sum of the own vortex field 
and the field of image vortices, situated inside the cav- 
ity. Also, using the conformal transformation technique, 
pinning potentials for more tricky columnar cavities have 
been derivedi&£ However, in the calculation of the pin- 
ning potential for non-circular defects only the field of 
the image vortices has been transformed, while the mod- 
ification of the own field of the real vortex has not been 
taken into account. 

In order to analyze small defects, a more complex ap- 
proach is required. For temperatures close to the su- 
perconducting critical temperature the Ginzburg-Landau 
(GL) approximation is a natural choice. It has been ap- 
plied for the analysis of vortex pinning on columnar de- 
fects of different shape and nature. Blatter at al4 es- 
timated the pinning potential for a small defect using 
a simple variational procedure, taking into account only 
the suppression of the order parameter inside the defect. 



Maurer et alM- calculated numerically the pinning energy 
for a vortex centered on a circular insulating or metallic 
inclusion. In Ref. [13 vortex interaction with a cylindri- 
cal hole with a square cross-section has been analyzed 
numerically. In Ref. the critical current for a vortex 
lattice pinned on a set of defects with reduced critical 
temperature has been determined using a variational pro- 
cedure and numerical simulations. Yet, exact analytical 
solutions of the GL equation have been lacking so far. 

In this paper, within the GL theory we consider the in- 
teraction of a vortex with a small cylindrical cavity or in- 
sulating inclusion with the characteristic size of the cross- 
section £ D <C £, where £o is the zero-temperature 
coherence length (in the case D <C £o the correct descrip- 
tion can be obtained only on the basis of a microscopic 
theory- 1 -"^). We present the exact pinning potentials in 
terms of the unperturbed vortex order parameter for a 
circular and elliptic defect. For the circular cavity the 
depinning current is also determined. For the treatment 
of a small defect with arbirary cross-section we propose 
an electrostatic analogy, which serves as a counterpart to 
the mentioned above analogy between the London the- 
ory (valid for large defects, D ^> £) and electrostatics. 
Finally, using the conformal transformation method de- 
veloped in Ref. |7j, we derive the pinning potential for an 
elliptic cavity within the London approximation, taking 
into account the modification of both the own vortex field 
and the image field. 



II. VORTEX PINNING WITHIN THE 
GINZBURG-LANDAU THEORY 

A. Basic equations 

The starting point for our analysis is the GL equation 
for the order parameter ip — \ip\ e : 



e 2 VV-V' + n " 1 |V'| 2 ^ = ! 



(1) 



where no is the concentration of Cooper pairs in the bulk. 
The GL parameter is assumed to be large, so the vector 



potential can be neglected^ Let us put the origin of co- 
ordinates inside the insulating defect, or cavity, and the 
z-axis along the vortex axis. Then if the vortex axis is 
parallel to the generatrix of the defect and perpendicu- 
lar to the transport current, the order parameter does 
not depend on z. This is the case that will be analyzed 
further. 

Equation (fTJ) is supplemented by two boundary condi- 
tions, specifying the normal derivative of ip at the insu- 
lating defect border and the transport current density j tr 
far from the vortex core: 
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Here S is the defect cross-section, dS denotes the border 
of S, n is the outward unit normal to dS, e is the electron 
charge, and m is the Cooper pair mass. 

Since our system must contain one vortex, an addi- 
tional condition for the order parameter phase arises: 



j> V6d\ = 2tt, 



(4) 



where integration is performed over a sufficiently large 
contour surrounding the defect. 

We expect that for some current jd Eqs. |T]) - (@| can 
be solved when j tr < jd, and no solution exists when 
jtr > jd- Then it is natural to consider jd as a depinning 
current. 

We will solve Eqs. (TJJ) - ^ in the case j tr < jd- If the 
transport current is much smaller than the depairing cur- 
rent, the order parameter has the following asymptotics 
at infinity: 



4, 



n e 
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where tp is the polar angle, and q = mj tr /2e?mo. This 
asymptotics can be derived from Eqs. ([TJ, ([3]) and (HJ 
if one expands tp in powers of p~ x and neglects terms 
proportional to j^. Now we make some assumptions 
concerning the behavior of the order parameter in the 
vicinity of the defect. 

(A) tp(p) reaches its asymptotic behavior at sufficiently 
small distances from the origin: tp s» y / 7^oe ^¥, • e 4qp 
when p > R, where R is some radius in the range 

(B) The vortex is weakly distorted by a small defect 
and a small current. This means that the solution 
of Eqs. ((T|) - (Q} can be presented in the form ip — 
tpo + ipi , where ipo corresponds to an unperturbed 
vortex shifted from the origin by a vector L (see 
Fig. [1]), and ip% is a small perturbation: IV'iCp)! •C 
^/n^ when p <C q~ x ■ 




Defect 



Vortex 



FIG. 1. The cross-section of the system: in the presence of a 
transport current jt r the vortex is shifted by a vector L with 
respect to the defect. 



The assumption (B) is justified by the fact that the un- 
perturbed vortex corresponds to a local minimum of the 
free energy, so large distortions are not energetically fa- 
vorable. Both statements (A) and (B) can be verified by 
numerical calculations. 

Let us write down the equations for the function ip\. 
If we linearize Eq. (TTJ we obtain 

- e 2 V 2 ^ - Vi + 2ng 1 |V> | 2 fa + no VoV>i - 0. (6) 
According to the statement (A), 

^i(p) = ip(p) - MP) » s/mie iv (<ip), R<p<^q- 1 . 

(7) 

The boundary condition at the defect border for tpi fol- 
lows from Eq. ©: 



(VV>i + V^o)n 



= 0. 



(8) 



as 



( 5 ) Thus, Eqs. © - © are to be solved. 



B. Variational derivation of the pinning potential 

In this section, for the reader's convenience, we present 
a relatively simple, but not rigorous derivation of the 
pinning potential. A more detailed and careful analysis 
is given in section HTCl 

We will determine the free energy F per unit length of 
a vortex shifted from the origin by a vector L: 
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where H c is the thermodynamic critical field satisfying 
the relation 



H 2 n 
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In the zero-order approximation F equals the free energy 
of an unperturbed vortex: 



F = 



HI 



e\VM 2 -\M 2 + n^-^)cP P . 



The pinning potential equals the difference between the 
exact free energy F and F a : U P (L) = F(L) — F . This 
difference consists of two terms - Aiq and AF2. The 
first term is connected with the suppression of the order 
parameter inside the defect: 



AFx 



-S 



Hi 



Anno \ ' ' ' TU ' ' 2 

\ / P =o 

Note that this expression is valid only for very small 
defects with the characteristic size of the cross-section 
D -C £. The second term is connected with the distor- 
tion of the order parameter outside the defect: 



AF 2 = 
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We substitute here ip = ipo + ipi, where ipi is a small 
perturbation satisfying Eqs. ([6]) and (jHJ in the vicinity 
of the defect and decaying to zero on a scale p ~ £ (if 
we use the function ipi satisfying Eq. © in the whole 
xy plane, the component AF2 will diverge to positive 
infinity). Close to the cavity the characteristic length 
scale for ipi is D <C £, hence, we can neglect terms of the 
order of l^il 2 as compared to the term £ 2 |VV>i| 2 : 

- Jos bK(VVon) + ^i(VV^n)] di) 

Using Eq. (jHJ and applying the Gauss theorem, we trans- 
form the right-hand side of the last relation as follows: 



AF? 



eat 



Jp$S 



(|VV>i| 2 + ip*i VVi + VW 2 Vi) d 2 P 



From Eq. © we find that 

|ViVVi| ~ |^i| 2 /e 2 « |VVi| 2 

hence 



(10) 



AFo 
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|VVi| 2 d 2 p- 



(11) 



p<ts 



An explicit expression for AF2 in terms of ?/>o will be 
given below (see Eq. (|25|)). 



C. Force balance equation 

In this section we derive the solvability condition for 
the system © - ©. Our derivation closely follows the 



computations from Ref. [lj which were used to determine 
the viscous drag force acting on a moving vortex. 

First, we introduce the auxiliary function ipd — dVipo, 
where d is an arbitrary constant unit vector, ipd satifies 
the equation 

- e 2 V 2 V^ - ipd + 2V 1 |Vo| 2 i>d + n^VoVS = 0. (12) 

Let us multiply © by ip*, and subtract Eq. (p~2|) multi- 
plied by ip* from it. When we add the complex conjugate 
to the resulting equation we obtain 

div (-#|VV>i + ViVV'd - ^VV* + Vi v^) = 0. 

We integrate the last relation over the region p ^ S, 
p < R and apply the Gauss theorem: 

S p =r (-^Vi/.i + ipiViPd - i><Ni>i + ipivipa) mde 

~ Ids (-^rf v ^i + V'iV^ - ipd^Vi + VWd) nd£ = 0, 

(13) 

where ni is the outward unit normal to the circle p = 
R. The first integral can be calculated with the help of 
Eq.0: 



f p =H (-V>3Wi + ^iV^ - + VWd) nd£ 



2i7m 
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(14) 



where Zo is the unit vector along the z axis. 

The second integral in Eq. (fT3| can be transformed 
using Eq. © and the Gauss theorem: 

Ids (~^ v ^i ~ i><N<ipl + ^iV^5 + </>i Wd) nd£ 



5 -div (V'5VV'o + V'dVV'5) 



p=0 



+ /as WiV^(O) + l/W d (0)) ndl (15) 

Here and further we neglect terms which are much 
smaller than no-D 2 /£ 3 . 

In order to proceed we have to determine the value of 
ipi at the defect boundary. It has been noted in Sec. Ill Bl 
that tpi ~ JuqD^ 1 . Hence, 



d 2 ^/dx 2 ~ ~ v^o^r 1 » \h r 2 . 

This implies that near the defect we can use the Laplace 
equation 



vVi = 0. 



(16) 



instead of Eq. © . The boundary condition can also be 
simplified: 



V^in 



= -V^ (0)n. 



(17) 



as 



Such simplification is acceptable since we are not inter- 
ested in small corrections of the order of ^fn^D 2 /£ 2 to ipi. 
Equations (fTB| and (| 1 T[) are equivalent to an electrostatic 
problem where ip\ plays the role of the electric potential 



4 



of a charged cylinder. Note that these equations can not 
be derived within the electrostatic approximation for the 
London theory^ where variations of the superconducting 
phase are taken into account, but the order parameter 
modulus is assumed to be constant. 
The relation 

I V?AincZ£ = 0, 

JdS 

can be interpreted as a vanishing total "charge" of the 
cylinder. It provides that a solution of Eqs. ([T6|) and 
(fTTf exists that decays like p^ 1 at infinity. This solution, 
which we denote as ip^ , represents the irregular part of 
-01 : it has singularities inside the defect. We define the 
regular component of ipi as 



= ipi- "01' 



(18) 



It is proved in the Appendix that the contribution of 
ipi to the integral in the right-hand side of Eq. (fT5|) is 
negligible. 

Combining Eqs. (fT3|) - ([15} and using the fact that tf>o 
satisfies Eq. (TTJ) we obtain 



^ (d [z ; j te ]) + S(dV)(e |V^ | 2 - |0 O | 2 + ^ 



p=0 



+ ? fas (VfW^O) + 4 d) *VM0)) ndi « 0. (19) 

Owing to the linearity of Eqs. (|T6|) and (Ti"7|) the solution 
can be presented in the form 



# = (g(p);V^(0)), 



(20) 



where g is a real vector field defined by the relations 



V 2 g = 0, 
Then 



(nV)g 



OS 



0. (21) 



J as (<V^(0) + < } *V^(0)) nd£ 
= (dV) (v0oGVV^ 



p=0 

where G is a real symmetric matrix with components 

G i:j = I g^dt = I V gi Vg 3 d 2 p. (22) 
Jas J p^s 



Equation (|19j) transforms into 



Here 



-^(d[z o ;j tr ])-(dV L )C/ p = 0. 
nhc/e is the flux quantum, Vl = d/dli, and 



4-7rno 



2;>„ 



p=0 



? 2 4^tV^o(O)GV0*(O). 



(23) 



Since d is an arbitrary vector, it can be dropped, and we 
finally obtain the force balance equation, connecting the 
vortex displacement L with the transport current j tr : 



ZOlJtr 

C 



V L C/ p = 0. 



(24) 



Here, the first term is the Lorentz force and the sec- 
ond term is the pinning force: F p — — VlC/ p . Thus, we 
may conclude that U p is the pinning potential. It can 
be proved that this definition of the pinning potential is 
identical to the one given in Sec. Ill B[ if ip\ in Eq. (jTTJ) 

is replaced by ip[ d ^ . Thus, 



AF 2 



V^ (0)GV^*(0). 



(25) 



Before we determine some pinning potentials explicitly, 
we would like to note that our consideration can be easily 
generalized for the anisotropic case. Indeed, generally, 
the GL free energy can be presented in the form 



F=-&- f U 2 

4-im J I 



dx 



t 2 



dz 



I 2 + m 



(26) 



where £ y and £ z are the coherence lengths for dif- 
ferent directions. The scaling transformation x = x, 
V = y£,x/£,y, and z = z£ x /£ z reduces the free energy to 
the isotropic form. Thus, we again arrive at Eqs. ([T]) - 
©■ 

Now we consider two types of defects. 



D. A circular defect 

Let the defect be a circular cylinder with the radius a. 
When the origin is placed on the axis of the cylinder, the 
decaying solution of Eqs. (fT6|) and (|17|) is 



n 



(rf) 



a 2 (V^»o;p) 



(27) 



and the pinning potential is 



Up(L) 



2„2 



2? |V0 o | 2 -|0o| 2 



4no \ 2no 

\ / p=0 

(28) 

The function rpo can be determined numerically from the 
GL equation ((6]). We are not going into details of these 
calculations here. A detailed numerical analysis of this 
function can be found in Ref. [l5l 

Equation (|28[) allows us to determine the pinning en- 
ergy, E p , and compare it with a numerical result from 
preceding paper. According to our calculations, 



E„ 



U p (oo) - U p {0) = 0A7H 2 i 



which coincides with the numerical value given in Ref. © 
up to a factor of the order of unity. 



5 



The profiles of the pinning potential and the pinning 
force are plotted in Fig. [5] The pinning force reaches 
its maximum at L — L CI — 0.84£, where F p — F cr = 
0.252H 2 a 2 /£_. When j tT > cF cr /|0 o |, the force balance 
equation (p4|) has no solutions, hence 



3d = cF a / \<fo\ = 0.252 



2„2 



(29) 



is the depinning current. When j < jd Eq. (1241) has 
two solutions due to the nonmonotonic behavior of the 
function F p (L), but the solution with the larger vortex 
displacement is thermodynamically unstable. Indeed, it 
can be easily proved that it corresponds to a saddle point 
of the correction AF to the free energy of a vortex con- 
nected with the presence of the defect and the transport 
current: 



AF = U v 



[(ZO X j tr )p] 



To sum up, deppining from a circular defect occurs 
as follows: as the transport current increases from zero 
to jd the vortex displacement with respect to the origin 
increases from zero to L cr ; as the current is increased 
further the vortex is carried away from the pinning cite. 



E. An elliptic defect 

Consider an elliptic defect with the cross-section 



x 2 y 2 



b 2 



< 1, 



where a > b. We shall determine the vector held g and 
the pinning potential. It is convenient to use the elliptic 
coordinates ((,T]): 



= V 'a 2 — b 2 cosh £ cos rj, 



V 



V a 2 — b 2 sinh £ sin 



The border of the defect corresponds to the value C = Co, 
where 



sinh Co = 



cosh Co 



v a 2 -b 2 Va 2 - b 2 

Eqs. ([2"Tj) in the new coordinates read 

= 0, (30) 



d( 2 " r drj 2 ~ u ' 



g 



C=Co 
The solution is 



= —b cos rj, 



= —a sin 77. (31) 



C=Co 



9x 



be Co ~ c 



cos rj, 



g y = ae^° ''sin 77. (32) 



Using Eq. (l22l , we obtain the components of the matrix 
G: 



Tib , Gyy — ira , G xy — G y 



0. 



(33) 




FIG. 2. Profiles of the pinning potential (a) and the pinning 
force (b) for a circular defect. 



Hence, according to Eq. (|23[) . the pinning potential is 
a6(e 2 |V^ | 2 -|^o| 2 + ^ 



u =-SL 



+fb 2 



dx 



+ea 2 



dy 



(34) 



p=0 



The potential well for the vortex now does not have cylin- 
drical symmetry. As a result, the vortex displacement L 
and the depinning threshold jd depend on the direction 
of the transport current. 



III. COMPARISON WITH THE LONDON 
THEORY. 

In a preceding paper— the pinning potential in the pres- 
ence of a circular and elliptic cavity has been derived 
within the London theory. The applicability condition 
for this approach is that the defect size be much larger 
than the temperature-dependent coherence length, i.e., 
D » £(T), so it may seem that the results from Ref. 
|fj can not be compared with ours. However, the condi- 
tion D <C £ is not essential for our calculations. Indeed, 
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instead we may demand (i) |-0i| <C ^/uq and (ii): the 
quantity X?ipo should be approximately constant in the 
area occupied by the cavity. These two conditions are 
satisfied when 



D > £ and L > D, 



(35) 



so for a large defect and large vortex-defect distance our 
pinning potential should coincide with the one obtained 
within the London theory. 

For a circular defect our pinning potential (|28p at large 
vortex-defect distances, L>o, equals 



(36) 



This expression is in good agreement with Eq. (5) from 
Ref. [6j. For an elliptic hole Eq. (|34f yields in the L>s 
limit 



Up ~ {attXJ 2 I L 2 



1 / ab b 2 L 2 a 2 Ll 



L 4 L 4 
whereas the potential from Ref. is 

, \ 2 

a + b 



UttA J 



1 



-const, (37) 



(38) 



which, obviously, does not coincide with (|37|) . Below we 
will explain the reason of this discrepancy. 

The derivation of the interaction energy between a vor- 
tex and a cavity in the London approximation is based 
on the equation 



Uu 



0O^im(L) 

8tt 



(39) 



where h- ml is the z projection of the field created by image 
vortices. This field can be expressed as 



him h z 



2nX 2 



■In 



A 2 



C-Co 



(40) 



where £ = x + iy, Co = L x + iL y , and h z is the full 
magnetic field, satisfying the Poisson equation 



V 2 h z = -^S(p-L) 



(41) 



The second term in the right-hand side of Eq. (|4U|) rep- 
resents the own field of the vortex with the opposite sign. 
For a circular defect with the radius ao the image field 
at the position of the vortex is 



2nX 2 



In 1 



ICo 



(42) 



In order to obtain the magnetic field in the presence of 
a non-circular defect, we may apply a conformal trans- 
formation w — w(£) to the £ plane. Since the form of 



Poisson's equation is not modified by such a transforma- 
tion, the field distribution in the w plane is given by 

h z (w) = h c z (((w)), 

where h c z (£) is the solution of Eq. (|4"Tj) in the presence of 
a circular defect. Using the definition (|40[) of the image 
field, we obtain 



frim(wo) = hf m (((w )) 



2ir\ 2 



111 



A 



2-rrX 2 



hi 



, (43) 



where wq = u>(£o) specifies the position of the vortex in 
the w-plane. Hence, the pinning potential equals 



In 1 



R 2 



icK)r 



In 



d£ 
dw 



Here, the first logarithmic term originates from the trans- 
formation of the image field (1421) , while the second term 
is connected with the modification of the own field of 
the vortex. In Ref. la this term has not been taken into 
account. As a result, the isotropic potential (|55]l has 
been obtained. In order to determine the correct pinning 
potential for an elliptic cavity, we apply the modified 
Joukovski transformation: 6 



Uu 



\4tt\ J 



b £ 



ba 



(45) 



In 1 



a+b 



In 



1 



y/w 2 — a 2 +fr 2 



w-\-y/w 2 — a 2 -\-b 2 

const. 



(46) 



For \w\ S> a this expression coincides with our result, 
obtained within the GL theory (see Eq. (|3"?| ). 



IV. CONCLUSION 

By solving the Ginzburg-Landau equation, we devel- 
oped a method to determine analytically the pinning po- 
tential for a vortex interacting with a small cylindrical 
cavity. This method has been applied to a circular and 
elliptic defect. In the latter case, the pinning potential 
appeared to be anisotropic, as one would expect. Also, 
we recalculated the pinning potential for an elliptic cavity 
within the London theory, using the conformal transfor- 
mation technique 23 - considering the modification of the 
image field as well as the transformation of the own field 
of the vortex, which had not been previously taken into 
account. Our results obtained within the GL and Lon- 
don theories agree well with each other in the range of 
parameters, where both approaches are valid. 

All our previous consideration has been related to the 
case of a vortex strictly parallel to the defect. This as- 
sumption is obviously satisfied in thin superconducting 



7 



films, where the transport current is distributed almost 
uniformly over the film thickness. However, it has been 
claimed^ that in a bulk superconductor depinning is 
likely to occur due to vortex kink formation in a surface 
layer with the thickness equal to the London length. Still, 
if the vortex radius of curvature is large as compared to £, 
and the vortex axis makes a small angle with the defect 
axis, our approach should give reasonable estimates of 
the pinning energy and of the local vortex dispalcement 
with respect to the defect. Moreover, in the presence of 
sufficiently small transport current a bound state should 
occur, when the vortex core is outside the defect (at least 
in a surface layer - in bulk superconductors) . This bound 
state should be observable, for example, using scanning 
tunneling microscopy. 

Our results may be useful for estimations of the depin- 
ning current and for manipulations of the critical current 
anisotropy in superconducting materials. 
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Appendix A 



In this appendix we will demonstrate that the con- 
tribution from the function tpY (see Eq. (JT8J) ) to the 
pinning force is negligible. It is sufficient to prove that 



the absolute value of the integral 

I I (#VV>K0)+#*VV>d(0))nd£ 



as 



is much smaller than uqD 2 /^ 3 



The function tp\ has the following properties: 
*0, 



VV£° =0, p < r, (Al) 



as 



where r is a quantity of the order of the coherence length. 
Let us introduce an auxiliary function v defined by the 
relations 



V 2 v = 0, 



as 



The properties of this function are identical to those of 
ip[ . it is of the order of y/noD^ 1 at the defect border 
and decays like p~ x at infinity. For a smooth defect v ~ 
y/n^D 2 /^p. Now we make some simple calculations: 



a P 



dC. 



Since ip^' w tpi when p ~ r, 

i = r 1 



dv dipx 

Ipl- V^r— dl+CC, 

v dp dp 



(A2) 



where c.c. denotes the complex conjugate. According to 
statement (B) from section llTAl when p = r , <C y/no 
and \dipi/dp\ <C \Zno/£, since the characteristic length 
scale is £. Then it follows immediately from Eq. (|A2[) 
that \I\ « n D 2 /£ 3 . 
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